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We demonstrate by an explicit model calculation that the decay of entanglement of two two-state 
systems (two qubits) is governed by the product of the factors that measure the degree of decoherence 
of each of the qubits, subject to independent sources of quantum noise. This demonstrates an 
important physical property that separated open quantum systems can evolve quantum mechanically 
on time scales larger than the times for which they remain entangled. 
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Entanglement of quantum-mechanical states, referring 
to the nonlocal quantum correlations between subsys- 
tems, is one of the key resources in the field of quan- 
tum information science. Many protocols in quantum 
communication and quantum computation are based on 
entangled states When one considers practical ap- 
plications of entanglement, the coupling of the quantum 
system and its subsystems to the environment, resulting 
in decoherence, should be taken into account. It is known 
^, 3J that entanglement cannot be restored by local op- 
erations and classical communications once it has been 
lost, so understanding of the dynamics of decoherence of 
entanglement is of importance in many applications. 

There are two basic issues in the physics of the loss of 
entanglement by decoherence, that, while intuitively sug- 
gestive, thus far have allowed little quantitative, model- 
based understanding. To define them, let us refer to 
two subsystems, S^^^ and 5^^^ , of the combined system, 
S. The first property of interest is the expectation that 
when the systems are separated in that they are subject 
to independent sources of noise, e.g., when they are spa- 
tially far apart, then the decoherence of entanglement is 
faster |j, |M |y{ than the loss of coherence in the quantum- 
mechanical behavior of each of the subsystems. Thus, 
the subsystems can for some time still behave approxi- 
mately in a coherent quantum-mechanical manner, but 
without correlation with each other. 

In order to define the second property of interest, let 
us point out that the definition of "decoherence" of an 
open quantum system is not unique. One has to consider 
the overall time-dependent behavior of the reduced den- 
sity matrix of the system, obtained for a model of the 
environmental modes, which are the source of noise and 
are traced over. This time dependence can involve an os- 
cillatory behavior corresponding to the initial regime of 
approximately coherent evolution, with frequencies de- 
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termined by the energy gaps of the system (which can 
be shifted by the noise). At the same time, there will 
be irreversible, decay- type time dependencies manifest 
for larger time scales, which can in many cases be iden- 
tified with processes such as relaxation, thermalization, 
pure decoherence, etc^ that represent irreversible noise- 
induced behaviors 0, i, lE, 10, H .12, 1^ 14, 15., 16] . 

One, by no means unique, way to quantify the degree 
of loss of coherence is by the decay of the absolute values 
of off-diagonal elements of the reduced density matrix. 
This definition is only meaningful at relatively late stages 
of the dynamics, when the density matrix has already 
become nearly diagonal in a basis favored by external and 
internal interactions, and by environmental influences, 
e.g., for thermalization, the energy basis. More careful 
definitions of measures of decoherence are possible p^ . 
but we will use the off-diagonal-element nomenclature 
for clarity. Recent experimental NMR studies |0| have 
considered various "orders of coherence" that involve off- 
diagonal elements, for systems of up to 650 spins. 

The second property of interest, is formulated in this 
language as follows. For noninteracting and noncntan- 
gled subsystems, the density matrix of the whole system 
will be a direct product of the subsystem density ma- 
trixes. In this simple case, there will be far-off-diagonal 
density matrix elements of the system that will decay 
by a factor that is a product of the decay factors of the 
subsystem off-diagonal elements. Specifically, if the large 
time decay is exponential, then the decay rates will be 
additive .19. .20.]. 

A related "additivity" property has been mathemat- 
ically explored for certain measures of initial decoher- 
ence (17i| . for entangled subsystems. Recently, explo- 
ration of the following physically very suggestive ques- 
tion has been initiated |5(: If we know the suppression 
factors, < S'^^'^'^ < 1, that roughly measure decoher- 
ence for the two subsystems, then are there any phys- 
ically meaningful quantities that are suppressed by the 
product 6'^^^ 6^^^ 7 The other suggestive alternative is that 
the "worst case scenario" for physically relevant loss-of- 
coherence measures of the combined system is suppres- 
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sion by the factor of mm{S^-^\ S^^''). The two alternatives 
are, of course, only approximate, qualitative statements, 
possibly for upper bounds for oscillatory quantities, be- 
cause we have not specified the precise measures to use, 
nor the dependence on (or maximization of the decay rate 
over) the initial conditions. 

In this work, we show by an explicit calculation for a 
solvable pure-decoherence model of two qubits (two-state 
systems, spins-1/2) interacting with a bath of bosonic 
modes, that the measure of entanglement introduced in 
PH . is indeed suppressed by the factor J^^^J^^^. We fo- 
cus on the two-qubit system, because it is only for this 
simplest case that an explicit expression for a measure of 
entanglement called concurrence was obtained [U. Our 
study expands the recent works that considered sim- 
ilar properties for different models. We are able to derive 
explicitly the product of suppression factors result. 

For brevity, from now on we will use subscripts or 
superscripts r = 1,2 to label the spins (two- level sub- 
systems), Hg = A^cr^. Each spin interacts with a 

bosonic bath of modes Hp = X^fe'^fc^fe^^L which has 
been widely used @, 0, 0| as a model of quantum 
noise (we set h — 1). The interaction between the quan- 
tum systems and the environment is taken in the form 

= '^z X]fc {91*^1 + 5fe^fc^) • This choice, correspond- 
ing to [H'^jHj] = 0, leads to a solvable model and has 
been identified as an appropriate description of pure de- 
coherence 

We assume that there is no interaction between the 
qubits, so that the Hamiltonian of the whole system has 
the form H = J2r i^s + + ^i)- The main reason for 
this assumption is, of course, to have a solvable model. 
In addition, we point out that qubit-qubit interactions, 
either direct or those induced by the bath modes, can 
decrease or increase their entanglement. For the latter 
reason, we also assumed that the noise is uncorrelated 
at the two qubit locations, namely the bath modes are 
independent for each qubit (the most natural situation is 
when the qubits are spatially separated). 

The initial state of the two qubits, described by the 
density matrix ps(0), can be entangled. However, we 
assume [l^oJJJ that the qubits are initially not entangled 
with the bath modes. The overall initial density matrix 
is then 



p(0) = P5(0)®p^(0)®p|(0) 



(1) 



The reservoirs are in thermal equilibrium at the temper- 
ature T (with P = l/knT), 



The total Hamiltonian is time-independent, so the re- 
duced density matrix of the two qubits at time t > is 



Psit)^TrB [C/p(0)C/t] , ^ (3) 

where the evolution operator factorizes, U = e~*^* — 
UiU2- The trace over the bosonic modes of the two baths, 
Ttb in lO, can then be evaluated exactly by using the 
techniques of 0, . 

It is convenient to write the density operator psit) in 
the matrix form, 



{t)^{jh!\ps{t)\jhi) 



(4) 



where 7^ = ±1 has two indexes: r labels the qubit, while 
q simply indicates whether it marks row or column matrix 
element positions. The values -1-1 and —1 correspond to 
the spin states t and J,, respectively. 

After several straightforward transformations, ||3J) is re- 
duced to 



(5) 



where the coefficients are 

here AT is defined by iJJ = (J^Hj. Utilizing the identities 
from |7lll6l|. we find an explicit expression 

T^rT^ =exp[-G.(i)(7r-72^)'] , (7) 
where Gr (t) is the well-studied spectral function , 



Gr(t) = 2 



^sin^^coth^. (8) 



PB(0) = n(l-^"''"^)^ 



(2) 



A general property of the pure-decoherence models 0| 
is that the diagonal elements of the density matrix will 
stay unchanged during the evolution. 

Utilizing the new variables = e^*'^'^* and = 
g-4Gr(i) density matrix can be written explicitly. 
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PS {t) = 



/ pV-'^ (0) 

^292^^'" (0) 

(0) 

VpigiP2g2pi^'"(0) 



P2924^'^^ (0) 

pI''^' (0) 

PiqiP*2q2pi^'^^ (0) 
PigiP^^'T^ (0) 



Ti.JT * TJ.,iJ, ■ ' 



pIIiPs 
P*iqiP2q2pl''-" (0) 

P2q2pii^'" (0) 



Pi'Zipij'''' (0) 
Plq2pi^-^^ (0) 



(9) 
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FIG. 1: Eigenvalues fii and /i2 as functions of rj, for several 
values of ^, with the prefactor |ap/(l + suppressed. 



To analyze the effect of decoherence on the entangled 
qubit states we use a measure of entanglement. The 
entanglement of formation [221 was historically the first 
widely accepted measure of entanglement. For a mixed 
state PS, the evaluation of this measure is related to min- 
imization over all the possible pure-state decompositions 
of PS, and even for a two-qubit system getting analytical 
results for this measure is a complicated problem. The 
concurrence \2l\ is a quantity monotonically related to 
the entanglement of formation, hence it may be used as 
a convenient substitute for it. Given a pure or mixed 
state, ps, of two qubits, we define the spin-flipped state 



PS = (o-y (Xi cTy) pg {ay (g) ay) 



(10) 



and the Hermitian matrix R{ps) = \/ \fPs PS\fPS with 
eigenvalues A^^i 234- Then the concurrence |21j is given 
by 
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C(/9s(0) = max{0, 2maxA, -^A^} . (11) 



Since we know ps{t) explicitly the evaluation of pi|) 
is reduced to finding the eigenvalues of a 4 x 4 matrix. 

For illustration, we considered the system of two 
qubits in a pure state which at time t = Q is = 

(Iti) + a I jt)) /yl + H^- Here the (complex) param- 
eter a characterizes the degree of entanglement. Under 
the influence of the quantum noise the system evolves 
from the state ps (0) — |V') (V'l to the mixed state 



PS {t) 



l + \a\ 



/ 

1 

piqiplq2a 
V 




P*iqiP2q2a* 



\a\' 




(12) 



To evaluate the measure of entanglement at times 
t > 0, we have to find the eigenvalues of the matrix 
R, which can be obtained from the eigenvalues of the 
product Ps (t) Ps (t)- The latter eigenvalues are 
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1+^ cos2r/±2^cos77y 1 — ^2 sin 77 



(13) 



and /.t3_4 = 0. Here ^ = e-4[Gi(t)+G2(t)] and r/ = 
2 (^A^ — A^) t. Then the eigenvalues Xi = y/pi , and as a 
result the concurrence takes the form, 



Cips {t)) 



/pi 



(14) 



The eigenvalues pi^2 are shown in Fig. 1. For example, 
for a simple case of identical qubits, A^ = A^ , we have 
77 = and Xia = |a|(C ± 1)/(1 + |aP)- As a result the 
concurrence is 



2\a\ 



-4[Gi(t)-(-G2(t)] 



(15) 



This establishes the product of the suppression factors 
property alluded to in the introduction, because it is 
known Till that each of the exponential factors e"^'^!'^'^*^ 
measures the decay of the off-diagonal matrix elements 
when each qubit is isolated from the other, but exposed 
to its own bath. When the qubits are not identical, one 
can prove that for any t > 0, 



(16) 



so that the product of the factors property applies as an 
upper bound. The recent Markovian-approximation re- 
sults [^1^, appropriate for large times, have yielded an 
interesting observation that for some initial conditions 
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the concurrence, unlike coherence, can drop to zero in fi- 
nite time (231, (2^ 125(1. We have not explored this property 
within the pure-decoherence scheme considered here. 

In summary, we connected two important issues in the 
studies of entanglement and dccoherence, namely, for a 
solvable pure-decoherence model, we confirmed that the 
decay of entanglement is approximately governed by the 
product of the suppression factors describing dccoher- 
ence of the subsystems, provided that they are subject 
to uncorrelated sources of noise. Our results also suggest 
avenues for future work. Specifically, for multiqubit sys- 
tems, one might speculate that similar arguments could 



apply "by induction." However, understanding of entan- 
glement is far from intuitive, especially when one con- 
siders more than two two-state systems. Therefore, for 
any definitive progress, one has first to develop appro- 
priate quantitative measures of entanglement, and try to 
quantify entanglement and dccoherence in a unified way. 
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